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Abstract 

In this paper we derive the scaling fields in c = — 2 conformal field theory associated 
with weakly allowed clusters in abelian sandpile model and show a direct relation between 
the two models. 
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1 Introduction 

There exists some phenomena which naturally show power law behavior, that is, without fine 
tuning any parameter, the system shows behavior similar to the critical point, in contrast 
with the usual critical phenomena, where you should fine tune an external parameter like 
temperature to arrive at the critical point. These kind of phenomena are said to have self 
organized criticality [Q. Sandpile jTJ |^] , surface growth |3| and river networks 0] are a few 
examples of such phenomena. 

The concept of SOC was first introduced by Bak, Tang and Wiesenfeld |X]. While many 
other models have been found after that, still the abelian sandpile model (ASM) is one of the 
simplest and most studied models. Despite its simplicity, it shows all the features that a self 
organized critical phenomena ought to present, so huge amount of work has been done on this 
model [5-18]. Many analytic results have been derived, for example the probability of different 
heights and many specific clusters are calculated explicitly El • Also the relation of this 
model to other known models has been discussed. First of all, there is a connection between 
ASM and spanning trees [7]. Then other models such as dense polymers, Scheidegger's model 
of river networks and g — s- limit of g-state Potts models arise |19|. 
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These statistical models, show conformal symmetry while they are at their critical point. 
So it is natural to look for a conformal field theory (CFT) which corresponds to these sys- 
tems. The CFT associated with these models is suggested to be c = —2, which belongs 
to a specific group of CFT's, known as logarithmic conformal field theories (LCFT's). In 
LCFT's, where correlation functions may have logarithmic terms in contrast with the or- 
dinary CFT's, there exist pairs of fields with the same conformal weight, which mix under 
conformal transformations 20} I21j. 

Maheiu and Ruelle jlOj have found a way to relate ASM to c = —2 theory. They have 
found some operators in the LCFT model which correspond to different clusters in ASM. But 
the correspondence is shown only through correlation function and it is not clear why one 
should take these operator. In this paper we'll address this question and find a direct way to 
derive the operators from the action of c = — 2, and hence connect ASM to c = — 2 directly. 
The paper is organized as follows: in section 2 we will briefly introduce the model and some 
analytical results including the result obtained by [TQ]. In the 3rd section we will derive the 
previously mentioned operators directly from the action of c = —2 model. 

2 Abelian Sandpile Model 

The Abelian Sandpile Model is defined on a square lattice of the size L x M. To each site i, 
a height variable hi is assigned. Height of sand at each site can take one of the values form 
the set {1,2,3,4}. So the total number of different allowed configuration is equal to 4 LM . 
The dynamic of the system is defined as follows: at each step a random site i is selected and 
a grain of sand is added to that site. If the new height of sand becomes more than four, 
the column of sand is called to be unstable and topples, that is, four grains will leave the 
site and each of them will be added to one of the neighbors. So the total number of sands 
is conserved during the toppling process except at the boundaries where one or more grains 
leave the system. 

The toppling process can be stated in another way, which will be more appropriate. If the 
site i becomes unstable, hj will be decreased by amount of Ay, that is hj — > hj — Ay where 



The matrix Ay is called toppling matrix. 

After a while the system reaches a steady state, in which it shows SOC. It has been 
shown that in this state, the number of different configurations the system accepts is det A 
and the probability of all of them are the same. These configurations are named recurrent 
configurations in contrast with the transient ones which can only appear in the first steps of 
evolution, where the system has not yet reached the steady state. Note that the number of 
recurrent configurations is fairly smaller than the total number of configuration as det A ~ 
3.2 LM . Determining whether a given configuration is recurrent or transient is a relatively 
hard question, though there exist tests to answer this question. The first observation is 
that some forbidden subconfigurations exist, that is if in a given configuration you find one 
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Figure 1: Modification of the sandpile lattice in the two simplest Clusters Sq and Si 

of these sub configurations then it could not be recurrent. The simplest example of these 
subconfigurations is two neighboring height one sites. 

One of the most important analytic results derived so far, is the probability of finding 
some specific clusters in the system while it is in the steady state. Dhar and Majumdar 6 
have found a subtle way to calculate such probabilities for the clusters named weakly allowed 
clusters (WAC's). WAC's are those clusters which by decreasing the height of any of its 
sites by one, it becomes a forbidden sub configuration. The simplest one is a single site with 
height equal to one. They have proved that the number of recurrent configuration with a 
particular WAC is equal to the number of recurrent configurations in another sandpile model 
with a modified toppling rule. In fact there are several different ways to modify the toppling 
rule. One of them is to remove all the connections of the WAC to the rest of the system but 
one as shown in figure 1 for two simplest clusters known as So and S\ respectively 1 . As the 
grains of sand are not allowed to flow through the disconnected bonds, both the condition for 
instability and the toppling matrix should be modified. The new toppling matrix is given by 
A' = A + B, where for each disconnected bond running from site i to site j, Bij = Bji = 1 
and if n bonds have been cut from site i we should set Ba = —n. As it is clear, the matrix 
B is nonzero only in the vicinity of the cluster we are dealing with. So, B is somehow local. 

There are some other ways to modify the lattice and the toppling matrix. One of them 
which is more appropriate for our future uses, is to disconnect the cluster completely from 
other sites of lattices. This scheme is better because it preserves the symmetries of the 
clusters manifestly. The rules for modification of the toppling matrix is more or less the same 
as before, the only thing you should take into account is that Ba should be set to —3 for the 
sites in the cluster. 

Using these methods, the probability of finding different WAC's has been found analyti- 



1 If you are considering several disconnected pieces of WAC's, you should do the same for each piece. 
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cally |3 EI]. To do this, one should compute both det A and det A'. Then the probability of 
finding the cluster S would be 

P(5) = ^^ = d6t(1 + GjB5) ' (2) 

where G = A -1 , and A s and B$ are the matrices associated with the cluster S. As an 
example, the probability of finding a site with height one (Sq) is found to be P (1) = 2(1 — 
2 /it) /it 2 ~ 0.074. These analytic results have been confirmed by different simulations |22 [ 1231 
I24j though it seems there are some small disagreements |25| . Additionally, by takeing two 
clusters of this kind, placing far away from each other, one is able to compute correlations 
between these clusters. Taking two clusters to be So, separating with the distance r, one 
finds that the dominant term in correlation functions at large distances, is proportional to 
r -4 and is independent of angular distance of the two sites. For other clusters, the radial 
dependance is again like r -4 , but you'll have angular dependence. 

On the other hand, one is interested in the continuum limit of the problem, that is, finding 
a field theory that describes the same model. As ASM is shown to be equivalent to spanning 
trees and hence to q — > of q-state Potts model, the most appropriate filed theory seems 
to be c = —2 conformal field theory. Different features of the c = —2 model is described in 
many articles [2E|. It has the simple Gaussian action 

s = - [ 8989, (3) 

7T J 

where 9 and 9 are complex grassman variables. Comparing correlation functions, Mahieu and 
Ruelle |10j have shown that one can assign a suitable operator in c = —2 theory to each WAC 
in sandpile model. A table of operators versus clusters for some simple WAC's can be found 
in ^Hj. Actually, they have extended ASM to a model with dissipation in the following way: 
the threshold beyond which the column of sand becomes unstable is increased to x (x > 4) 
and during toppling process x sands are removed from the unstable site, but only four of 
them is added to the neighbors. So, x — 4 grains of sand are dissipated during each toppling. 
As a consequence, the action of c = —2 theory should be modified to 

s = ~ J (doBe + ^ee^j , (4) 

where m 2 is equal to x — 4. 

The operators assigned to different WAC's are very similar, in fact all of them are of the 
form 

(f> s {z) = - | A -.3656 + 5986: +B X : 8986 + 8089: +iB 2 -.8989 - 8989: +CP(S)^- 

(5) 

with A, B\, I?2 and C to be determined by comparing with the results of Dhar and Majumdar 
method. Also P(S) is the probability of finding the cluster S. The parameter C is present 
only in the massive theory and due to Mahieu and Ruelle it is very striking that C takes only 
integer values and is equal to number of sites of the cluster S. 
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In the next section, we will show a more direct way to find these operators. Also the 
connection between ASM and c = —2 model will be established more precisely. 



3 Connection Between ASM and c = — 2 Theory 

The first step to connect the two theories could be considering their action and partition 
functions. The c = —2 action is given by equation (J5J), but in the case of ASM we have not 
yet introduce the action. This can be done by considering that in the steady state there are 
det A different configurations with the same probability. So, the partition function of the 
model is just equal to det A. One can write the determinant of an arbitrary matrix in terms 
of Gaussian integration on grassman variables, namely: 

det A = J ]J dOi Y[ d9j exp (j^ Oi^ijGj) • (6) 

As the matrix A is the discrete version of Laplacian operator, it is clear that the above action 
leads to the action © in continuum limit. This was noted before by Ivashkevich |27j . 

Now we turn on to scaling fields assigned to WAC's. The starting point is equation (|2"|). 
where the probability of finding the cluster S is given. Again, the determinant of A can be 
written in terms of Gaussian integrals, so the probability P(S) turns out to be 



= detAg = /rfflirffljexp(Efl«A..flj) ^ j de . S . exp ( E (fl. A ^ + 9^9,)) 
[ ' det A Jd9 i d9 j exp (J29 i A ij 9 j ) J d9 { d9j exp (£ 0^0 j) ' [ ' 

Looking carefully at the above relation, one is able to derive the proper discrete field for the 
cluster S: 

<ps = exp 

The next step is to translate this field to continuum language. This translation should be 
done very carefully. Using the definition of the matrix B, the continuum version of the term 
inside the exponential is calculated easily. For the clusters shown in figure 1, we have 

X>£<i%ls=s oc 8989 + 8989, 

J2 e i B ii e ~j\s=s 1 oc 3(8989 + 8989) - (8989 + 8989) + .... (9) 

The exponential can be expanded in powers of 9 and 9, but since these variables are grassmans, 
the series ends at quadratic term. So, the operators we are interested in, are just the above 
fields, apart from an unimportant unity term in the series. 

The operator derived here for So is in complete agreement with the one derived in jlUj . 
you should just put the probability of finding height one to complete the proportionality. 
The field obtained for the cluster Si, however, is a little bit different, but not quite far away. 
The corresponding operator in JOJ is 



p Sl - 0.020143(0030 + 8989) - 0.006190(5050 + 8989), (10) 
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and the ratio of the two coefficients is about 3.25. Finding other operators corresponding to 
other clusters, one observes the same deviations, that is, although the ratios between different 
coefficient are not just the same as the one derived by jlUj . but they are very close to them. 
This shows that we are on the right path. 

To derive better result, one should do the process of continuation with more care. This 
means that we should first expand the exponential in terms of 9 and 9 before going from 
discrete to continuum. Let's consider a two point correlation function. Using Dhar and 
Majumdar method, the height correlation of two WAC's have the following form 

/ MidSj exp (E^iA^ + 9 i B° ij 9 j + 9^.9,)) 
Jd9 t d9 j exp (EOiAijOj) 

(exp (eiBffi) exp (fl^f/,) ) . (11) 

Expansion of the exponentials leads to several terms whose expectation values can be obtained 
using Wick theorem. We categorize the contractions in the following way: firstly, one may 
contract all the grassman variables in each of the clusters with themselves; secondly, only 
two connections are established between the two clusters 2 and so on. With every connection 
between the two WAC's, we will have a long range Green function. This means that the 
terms with more long range connections fall off more rapidly. As we are interested in the 
leading terms of the correlations function, only the first two terms could be considered 

P(S°,Sn = ((l + 9 i B? j 9 j + ...)(l + 9 k Bl l 9 l + ...)) 

= P(S°)P(Sn + (((1 + BiB^Oi + + 9 k B p kl 9t + ...))) , (12) 

where ((•••)) indicates that only the contractions with two long connections should be consid- 
ered. The first term, which comes from inter-contractions only, is simply equal to product of 
probabilities of the two clusters. However, in the second term one finds nontrivial correlations, 
and hence could be considered to contain the scaling operators we are looking for. 

The process of derivation of the operator associated with a given cluster can be summa- 
rized as follows. First we should expand the exponential exp(9iBij9j). Then contract all the 
variables in any possible way, leaving one 9 and one 9 uncontr acted. As an example if we 
have a term in the exponential of the form 9999 we should transform it to 

6>i6*2#3#4 — * ^1^2^34 — 9\9aG^2 + 9^9^G\2 — 6362G14. (13) 

Here, Gij is the Green function between the two sites i and j. The last step is to state the 
obtained expression in terms of 9, 9 at the origin and their derivatives. 

Doing all these together one is able to derive the desired scaling field. Though the pro- 
cedure is straitforward, it is very long and cumbersome. But fortunately, the method of 
expansion has a nice graphic representation as follows. First we put a point for every site 
involved in the matrix B. The contraction between 9i and 9j is represented with an arrow 

2 The number of connections should be even, because if you take it to be odd, odd number of 9 and 6 remain 
in each cluster and so the contribution of such term is zero 
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Figure 2: Two examples of many possible graphs related to the cluster So, and their corre- 
sponding expressions 

pointing from i to j, if i = j then we draw a loop. For the two sites whose variables are not 
contracted, one draws a dashed line from one to the other. Then for every line pointing from 
i to j, one puts the corresponding Green function, Gy multiplied by Bij. For the dashed line, 
one simply puts B^. The last rule is to put a minus sign for every loop with even number 
of lines. For example for the cluster So two of the possible graphs and the corresponding 
expressions are shown in figure 2. 

This equivalency makes the job more tractable. We draw all of the possible graphs and 
write the equivalent terms in the series. Then, By expanding 9i and Oj in terms of 9, 9 and 
their derivatives, we arrive at the scaling field of the specific cluster. The result for the cluster 
Sq is 

^(z) = - 2(7r 7 2) -.deae + Beoe ■., (14) 

which, apart from a factor of it, is the same as the one derived in JOJ- The 1/tt factor, comes 
in the very same way as it appears in the action (J2J). 

The result can be confirmed by doing it in another way. If we contract the remaining 9i 
and 9j, in the series (to produce an additional Gij), then the series will be equal to P{S), or 
equivalently det(l + BG). So, the term proportional to G^ in the expression for det(l + BG) 
is simply the coefficient of 9iQj in the above series. This method is much faster, specially if 
done numerically. For other clusters like S±, the corresponding scaling field is derived by the 
latter method and the results are in agreement with 

There remain some points to be clarified. First, we should note that any correlation 
function in c = —2 theory is zero unless you have the zero mode 69 in the correlation. So 
all the correlations like (4>s°4>sp) should be computed as (4 > s°4 > sp®Q) ■ This has a physical 
interpretation. The theory could not be defined on the whole plane, as the boundaries play 
an important role in ASM. In fact we should put mass term on the boundaries so that the 
grains of sand have the opportunity to leave the system. But if we would like to have the 
whole plane, it would be enough to put a field 09 at infinity. 

The other point is about the factor C in equation ifKj). By the method we derive the 
scaling field, it is quite natural that this factor should be an integer and the integer is just 
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the number of sites in the cluster. Because in the massive theory, the matrix B is a little 
bit changed: Ba is set to 1 — x instead of -3. So, for each site of the cluster we will have an 
additional term proportional to m 2 98, which is multiplied by the probability of the cluster 
when the contraction of other variables is done. 
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